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We show that maximally helical hypermagnetic fields produced during pseudoscalar inflation can 
generate the observed baryon asymmetry of the universe via the B + L anomaly in the Standard 
Model. We find that most of the parameter space of pseudoscalar inflation that explains the cos¬ 
mological data leads to baryon overproduction, hence the models of natural inflation are severely 
constrained. We also point out a connection between the baryon number and topology of the relic 
magnetic fields. Both the magnitude and sign of magnetic helicity can be detected in future diffuse 
gamma ray data. This will be a smoking gun evidence for a link between inflation and the baryon 
asymmetry of the Universe. 

PACS numbers: 98.80.Cq, 12.15.-y. 


Introduction .—Baryon asymmetry of the Universe 
(BAU) remains one of the greatest puzzles of cosmol¬ 
ogy. As was pointed out by Sakharov, any model that 
explains BAU has to satisfy three conditions: (1) baryon 
number non-conservation, (2) C and CP violations, and 
(3) departure from thermal equilibrium [T] . In fact, many 
of the early Universe problems were resolved within the 
inflationary paradigm. In this regard, one may wonder if 
BAU can also find its resolution within inflation. 

Since its proposal as a dynamical solution to the strong 
CP problem [5], axions provided a playground for rich 
phenomenology. These pseudo Nambu-Goldstone bosons 
(pseudoscalars) appear as a result of a broken global sym¬ 
metry. Besides, axion-like particles are abundant in string 
theory and are perfect candidates to build UV complete 
models of inflation, thanks to their radiatively stable 
properties. This stability is attributed to the existence 
of a flat direction that is protected by a shift symmetry. 
The symmetry eventually gets broken by nonperturba- 
tive effects, hence, the desired small values of the slow 
roll parameters arise in a technically natural way [3j. 

Pseudoscalars are naturally coupled to U(l) gauge 
fields through the dimension-5 operator (<f>/ f)Y ll „Y 111 ', 
which is inevitable from the effective field theory point 
of view. This term breaks the conformal invariance of the 
theory and leads to the production of cosmologically rel¬ 
evant fields^ [7] . The fields produced via this mechanism 
are coherent over the horizon scale and have maximal 
helicity [5]. Hence, they carry a non-zero Chern-Simons 
density, which breaks the macroscopic CP invariance of 
the Universe. 

In this letter, we point out that the change in the 
Chern-Simons density of the hypercharge field, h , pro¬ 
duced during inflation feeds into the baryon and lepton 
number anomaly in the Standard Model, which eventu¬ 


ally gets converted into the BAU. C invariance is broken 
since h enters the anomaly equation with opposite signs 
for different chiralities (see Eq. ([7j and Table [I]). In addi¬ 
tion, h provides the non-equilibrium condition required 
by Sakharov’s criteria and sustains the BAU. We find 
that natural inflation leads to overproduction of baryons 
for most of the parameter space compatible with the cos¬ 
mological data, hence this class of inflation models are 
severely constrained. In what follows, we use the metric 
= diag(l, —1, —1, —1) and set c = 1, h = 1, fcs = 1. 

Generation of the hypercharge field .—We consider a 
pseudoscalar inflaton $ coupled to the U(1) Y hyper¬ 
charge gauge field A M with held strength Y^ v : 

£ = m , (i) 


where V ($) is the inflaton potential, / is the axion con¬ 
stant, and a is the dimensionless coupling constant be¬ 
tween the axion and hypercharge held. The equation of 
motion of the gauge potential A M read^ 


/ d 1 2 „ 

2-V 2 -ayVx]4 = 0, 


( 2 ) 


where r is the conformal time and the prime denotes 
the derivative with respect to it.The generation of the 
hypercharge held can be envisaged by promoting the 
classical held A to an operator A, and then decompos¬ 
ing A into annihilation and creation operators aff as: 
A = Ea=± / (2 f) 3/2 [exAxa k x e lkx + h.c.], where the cir¬ 
cular helicity vectors e± obey the relations k ■ e+ = 0 
and k x e± = =pi|k| e ±- Substituting A into (J2J), we 
hnd that the mode functions A± satisfy the equation 
A’f + ( k 2 ± ak&/f) A ± = 0 (9j. 


1 See also Refs. mil for various mechanisms for generating pri¬ 

mordial magnetic fields. 


2 Here, we use the radiation gauge setting A 0 = 0 and V • A = 0. 
In this gauge, the electric and magnetic fields are E = —A'/a 2 
and B = VxA/a 2 , where a is the scale factor in the Friedman- 
Robertson-Walker Universe. 
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During the inflationary stage dQ/dt = <f>o is constant, 
where t is the cosmic time. The scale factor is given by 
a = —l/(Hi n { r), where f = a'(r)/a 2 (r) is the Hub¬ 
ble parameter during inflation. Hence, the equations of 


motion of the modes A+ read 


<TA± 


k 2 =F 2k ^ 

1 T 


= 


1 ± lcaAa dr 2 

0, where we defined the dimensionless parameter £ = 
a 9 /^° For A; |2£/r| the modes are in their vacuum. 
However, these modes will develop an instability when 
k ~ |£/r|. Depending on the sign of £, either A + or A_ 
modes will be amplified. For late times, \kr\ <C 2£, it is 
found that 0 


A+ = 


1 


a/2 k \ 2 £ a(r)H in { 


1/4 


e 7 r {-2 v / 2 {fc/[a(T)i4 in f] _ ( 3 ) 


Therefore, either A + or A_ will be amplified by a factor 
e 77 / and thus the produced hypercharge field is maxi¬ 
mally helical. This feature is an essential ingredient of 
the baryogenesis in pseudoscalar inflation. At this point, 
we emphasize that all the inodes produced during in¬ 
flation get diluted except the last mode that exits the 
horizon right before the end of inflation. This is also 
the mode that enters the horizon at the onset of re¬ 
heating and is the source for the BAU. In the above 
calculations, we neglected the backreaction of the gen¬ 
erated hypercharge field on the inflaton. This is a good 
approximation as long as the energy density stored in 
the generated field is less than the inflaton energy den¬ 
sity pinf = 3m 2 if? lf /(87 t). Using the mode decomposition 
described above [see also footnote <§] we find that the 
energy density stored in the hypercharge field is given by 
(E 2 + B 2 )- m { = 6!e 27r «tf 4 f /(2 19 7r 2 £ 3 ). Thus, we obtain 
an upper bound on the inflationary Hubble parameter 
H ini < 64^/tt/I 5£ 3/2 e _7r? m p . 

The hypermagnetic helicity is defined as T~L = J d 3 x A- 
B , which is proportional to the Chern-Simons num¬ 
ber. The rate of change of the helicity density is h' = 
—21imy_ >00 (l/U) f d 3 x E ■ B. As we will show below, 
the quantity of interest is (B ■ V x B). During inflation, 
the mode decomposition described above yields 


(B • V X B)inf = 


I 


d 3 fclfcl 3 
(27r) 3 


(l^ + | 2 


1 - 4 - 1 2 ) . ( 4 ) 


where the integral is over the comoving wave-vectors k. 
Using Eq. ([3]), cutting off the integral at k c = 2£ffi n f a(r) 
and setting one of the modes A± to zero, as only one of 
the two modes gets amplified, we fmcQ (B-Vx B) j n f = 
dp^H? n{ , where 1= 6.8 x 1CT 4 . 

As soon as reheating starts, the Universe becomes filled 
with a plasma of relativistic particles. The relevant mag- 
netohydrodynamics (MHD) equations are dtE = V xB 


3 Here, we choose the appropriate mode that results in a positive 
(B ■ V X B) i n f. This will ensure that we generate baryons rather 

than antibaryons. 


J,J = a(E+vx B), d t B = — V x E, where v is the fluid 
velocity of the plasma and er ~ 100T [lOj is its conductiv¬ 
ity. Neglecting the dtE term, which remains small in the 
MHD approximation, we obtain E = -'V x B — v x B , 
d t B = Vx(«xB) + -\7' 2 B. From these equations, we 
find [ll| E B = 4 B V x B, where we used BvxB = 0. 
On the other hand, in a plasma with a finite conductiv¬ 
ity, the evolution of hypermagnetic held as well as he¬ 
licity are governed by a competition between the dissi¬ 
pation, (l/cr)V 2 S, and advection, V x (v x B), terms. 
To this end, we recall the magnetic Reynolds number 
(the ratio between the advection and dissipation terms) 
lZ m = va/kp where k p is the physical wave-vector of the 
mode of interest, which is taken to be the last mode that 

exits the horizon at the end of inflation: k p = pA 

^ . .h 

where T r h is the reheating temperature. Assuming in¬ 
stant reheating, we have T r h = 0.25^/m/Hh/. Hence, 
we obtain lZ m = 400r>£ yfm p /Hi n {, which is much big¬ 
ger than unity for velocities v > y . In 

what follows, we assume that the cosmic fluid velocity 
is large enough that lZ m > 1 will be unavoidable un¬ 
til the electroweak phase transition. Hence, a turbulent 
flow will be generated, and the hypermagnetic field will 
no longer diffuse (see, e.g., Ref. [T2] for a review). Thus, 
h = -2 (E ■ B) = -(2 /a)(B ■ V x B) at a given epoch 
after inflation redshifts as: 


h 


2Xe 2 ^ (H ini y 
a£6 V a J ’ 


( 5 ) 


where we set aj n f = 1. In fact, it is this change in the 
helicity density that plays a pivotal role in baryogenesis 
by sourcing the B + L anomaly, as we show next. 

Chiral Anomaly in the Standard Model.— The Stan¬ 
dard Model fermions 



exhibit chiral anomaly; namely, the baryon and lepton 
numbers are anomalous in the Standard Model jlSj . The 
anomaly equation for a given fermion species / is: 


du^ = Ci^-Y UV Y» V + cl°pW“W a ^ 




167T 

C[^G b 


87 r 


87 r 


11V 


G b ^, 


( 7 ) 


where the coefficients C y , C w and C s are given in Table |TJ 
and a y , a w and a 3 are the hypercharge, weak and strong 
fine structure constants, respectively. 

In the vacuum, the rate of these anomalous 
baryon/lepton number violating processes are suppressed 
by e~ 2Sl) = e -47r / Q , where Se is the action of an instan- 
ton associated with this process, and a is the fine struc¬ 
ture constant of the corresponding gauge field. In partic¬ 
ular, the suppression factors are 0(e -164 ) and 0(e~ 55 ) 
for the weak and strong interactions, respectively. How¬ 
ever, in the early Universe, the dominant contribution to 
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Cy 

c w 

c s 

Q 

N c N w y Q 

Aft 

Aftv 

L 

Ny,yl 

1 

0 

Ur 

-N c y 2 UR 

0 

-1 

dR 

-NcVdn 

0 

-1 

Gr 

-yl R 

0 
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TABLE I. Coefficients C* in Eq. ([ 7 ]). The multiplicities 
Aft = 3 and Aft, = 2 take into account the color and weak 
isospin states of a given family of leptons and quarks, and the 
hypercharges are y Q = 1/3, y L = -1, y UR = 4/3, y dR = 
—2/3, y eR = —2. The charge conjugates Q c , L c , u R , d R and 
e c R have the same coefficients, C^, with all the signs flipped. 


the anomalous process is provided by sphalerons, whose 
rates can well exceed the expansion rate of the universe 
m- Therefore, the weak and strong sphalerons become 
efficient above the electroweak and QCD scales, respec¬ 
tively, and convert various fermions species into one an¬ 
other while preserving the gauge charges and the baryon 
minus lepton number, B — L. The effect of the hyper¬ 
charge sector U(1) Y is absent when the vacuum is triv¬ 
ial; namely, when both hyperelectric and hypermagnetic 
fields are zero, hence there is no change in the Chern- 
Simons number of this sector (see however Ref. [T5]). 
However, a change in the Abelian Chern-Simons num¬ 
ber, such as the creation of hypermagnetic field with 
net helicity, contributes to the anomalous processes. This 
effect is readily seen by integrating Eq. Q to obtain 
A N f = -Cf% f d 4 5 xE ■ B = ctjgAH, where AH is 
the change of the total helicity and ANf is the change in 
the baryon number. In fact, it is this integrated version of 
the anomaly equation that carries the physical informa¬ 
tion since the anomaly is related to the global structure 
of the theory. The integrated version tells us that AH is 
converted into ANf or vice versa. If the Universe during 
reheating were a bad conductor, then AH produced dur¬ 
ing inflation would survive the reheating to source ANf. 
If on the other hand, the Universe were a perfect conduc¬ 
tor, with infinite conductivity, then we would have E = 0 
and hence AH = A Nf = 0. This can also be seen from 


Eq. (|5j) as we find ft, = 0 by sending er —>• oo. However, 
the Universe has a finite conductivity, and thus, ft 7 ^ 0. 


In fact, ft is only suppressed by a factor ^ = 1 


1005 


H inf 

Dm D 


compared to the poor conductor case. 


In what follows, we study the effect of the maximally 
helical hypercharge fields produced during inflation on 
the evolution of the BAU taking into account the fi¬ 
nite conductivity of the Universe during reheating. To 
this end, we define ly as the asymmetry parameter for 

a given fermion / as 77 / = nf s ” /c ~ , where n/ is 

the number density, s = 27r 2 g»T 3 /45 is the entropy den¬ 
sity, 3 * = 106.75 is the number of the effective degrees 
of freedom in the primordial plasma for T > 10 3 GeV, 
and Hf is the chemical potential. We define the asym¬ 


metry parameters for each particle species in the Stan- 
dard Model as rj Q = ^N w N c fi Q T 2 , r) L = ±N W ^ L T 2 , 

Vur = gs^Vc, T]d R = 0^-/V c /i dH X' , T] eR = Q^fJ JeR T , 

and the baryon and lepton asymmetry parameters are 
ftbar = ^(Vq + t 1u r + Vd R ), V lep = Ng{V l + Ve R ), where 
N g = 3 is the number of generations. By taking an ensem¬ 
ble average of the anomaly equation ([ 7 ]) and taking the 
difference between the equations of the fields and their 
charge conjugates, we obtain the Boltzmann equation^ 
for ly (see for instance Ref. jT6] for the derivation of the 
Boltzmann equation for anomalous charges): 

^ = Fwfc+%) (8) 

- C£T b (t] q - rju R - Vd R ), 


where t is the cosmic time, and we used the relation 
= -4 (E ■ B) = 2ft. The rates T w ~ 25a 3 T [T7] 
and T s ~ 100 afT [T5] are the weak and strong sphaleron 
rates, respectively, and the coefficients C? are given in 
Table [I] In the radiation era we have H = 1/(2 1 ), or in 

terms of temperature T, H = jf -—, where D = ' 90 


Dm v ' > VVAA '“ av " ^ -y 87r 3 g* ’ 
In what follows, it will be more convenient to use the di¬ 
mensionless variable x = Dm p /T instead of the cosmic 
time t. After making the appropriate change of variables, 
the Boltzmann equation for the asymmetry parameter ?y 
of a given fermion can be written as a function of x: 


= ~ C yH ~ C'w7w(?7q + Vl) 

^ 7 s (fto - lu R ~Vd R ), 


(9) 


q- 2tt£ / tt \ 5/2 

where we defined 7y = , 7 w = 25a 3 , 

and % = 100a 3 . Having obtained the five Boltzmann 
equations for the five asymmetry parameters, we can 
readily integrate the system numerically starting from 
the reheating temperature that corresponds to some in¬ 
flation scale ffi n f all the way down to T = 10 TeV 
when the weak sphalerons shut of0 In Figs. |T] and |] 
we show the results for a range of the parameter space 
f° r 1 ^ £ < 5 taking zero initial values, r//(T r h) = 0. 
It is important at this point to emphasize that none of 
the above calculations are robust for £ < 1 , as in this 
case one can no longer trust the results for the generated 
hypercharge field [Eq. ©]- Th e plots in Figs. |T] and |]il- 
lustrate that the baryon asymmetry is over produced for 
most of the parameter space. Note also that one should 


4 The possibility that baryon asymmetry can be generated from 
helical primordial hypermagnetic fields via chiral anomaly in the 
Standard Model was pointed out in Ref. m- 

5 In fact, although the h term contributes to the evolution of rjf 
until the electroweak phase transition, T ew ~ 100 GeV, for sim¬ 
plicity we assumed that it shuts off at the same temperature 
as the weak sphalerons. We also neglected the Yukawa terms in 
Eq. d8l, which can only have a minor effect on the final rjf. 
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H inf (GeV) 

FIG. 1. Top panel r/ b ar vs £, the bottom panel y b ar vs Hi n f. 
For instance, the observed value of baryon asymmetry 7jbar = 
10~ 10 can be obtained for £ = 1 and Hj n f = 6.3 x 10*° GeV. 



FIG. 2. fl- inf vs £ for ? ?bar = lO' 10 . 


not extrapolate our results for values of ??bar > 1 as the 
perturbative calculation breaks down. 


Discussion .—A simple example of pseudoscalar infla¬ 
tion is natural inflation [3] in which a shift symmetry $ —> 
4> + C is broken down to a discrete subset <!>—>•$ + 27 t/ 
resulting in the potential V = A 4 [1 + cos(4>//)]. We re¬ 
call that the slow roll conditions are satisfied provided 
that to| < H ? nf . Using Friedmann’s equation 3fhpH? n{ = 
V — A 4 , where fh p = m v /y/&K, and m| = A 4 // 2 , we 
obtain / > \J3m p . The scalar fluctuation power spec¬ 


trum is given by = V 


1 + 7.5 x 10 ~ b V^ 


, where 


v 1/2 = df m- This power spectrum is probed by 
the CMB observations with amplitude given by the Cos¬ 
mic Background Explorer (COBE) normalization = 


25 x 10 _1Q . Using the equation of motion of the inflaton 
during the slow roll regime = dV/d<& = A 4 //, 

we find H m if — lCU 4 m 2 . Setting / ~ m p , we obtain 
H in f = 10 14 GeV. Such a large Hubble parameter dur¬ 
ing inflation will result in baryon asymmetry overpro¬ 
duction for values of £ > 1. One way out is that the 
coupling between the hypercharge gauge field and the in¬ 
flaton is very weak, a <C 1, such that no hypercharge 
field can be produced during inflation. Such fundamen¬ 
tally very small values of a appear to be contrived since 
one expects a > 1 as a consequence of the ’’gravity as 
the weakest force” conjecture m- Moreover, a consis¬ 
tent theory of quantum gravity disfavors values of the 
axion constant / > fh p . The latter problem can be solved 
within the framework of N-flation [2T]. In this scenario, 
one assumes that there are N different axions with con¬ 
stants ft = /single < iti v and that all these axions couple 
equally to the U(1) Y hypercharge gauge field such that 
effectively we have / = y/N /single > 0. Demanding 

that / < fhp one finds that N ~ 0(100) removes the con¬ 
flict with quantum gravity. On the other hand, lowering 
the inflation scale, H in f < 6.3 x 10 10 GeV, guarantees the 
production of the observed baryon asymmetry for £ > 1. 
Such low inflationary scales will require invoking curva- 
tons in order to respect the COBE normalization [ 221424 ] . 

One of the key predictions of natural inflation, apart 
from the cosmological data, is the maximally helical hy- 
permagnetic field generated via the dimension-5 opera¬ 
tor of the form (<&/f)Y pv Y ,lv , which is expected from 
effective field theory considerations. In this letter, we 
showed that the dramatic consequence of this coupling 
is the overproduction of baryon asymmetry that severely 
constrains models of natural inflation with large Hubble 
parameter, H ln f > 6.3 x 10 10 GeV. Note, however, that 
we assumed the primordial plasma to be turbulent be¬ 
tween T r h > T > 10 TeV. If the plasma ceases to be tur¬ 
bulent during the course of baryon number generation, 
the contribution from the hypercharge sector might be¬ 
come less efficient due to the diffusion of hypermagnetic 
field. Meanwhile, weak sphalerons can wash out the ex¬ 
cess baryon number. However, we note that this scenario 
is not likely as the magnetic Reynolds number tends to 
increase after reheating [T2|. To conclude, for parameters 
.ffinf < 6.3 x 10 10 GeV and £ ~ 1, the observed BAU can 
be achieved. Yet another prediction in this case would 
be the relic magnetic fields with right handed helicity 
(h > 0). It was pointed out in Refs. [I5}[25l|26] that the 
magnetic helicity is proportional to the baryon number. 
In this work, we explicitly showed this to be case. A re¬ 
cent analysis of the diffuse gamma ray data hints towards 
a global CP violation, which could be due to primor¬ 
dial magnetic fields with non zero helicity m- It would 
be a boon to find an observational correlation between 
the topology of these primordial magnetic fields and the 
baryon number. This will be a smoking gun evidence for 
a link between inflation and the BAU. 
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